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Abstract 

Macroscopic  plasma  diffusion  as  governed  by  the  simplest 
form  of  Ohm's  law  for  small  mean  free  path  Is  a  very  complex 
phenomenon  which  Is  by  no  means  described  by  the  classical 
"i/B  "  formula.   Diffusion  across  a  magnetic  field  is,  in 
general,  a  nonlocal,  nonlinear,  and  transient  combination  of 
skin  effect,  field  diffusion,  plasma  diffusion,  and  singular 
convection.   It  turns  out  that  several  orders  of  magnitude  are 
at  stake  in  balancing  these  competing  effects.   There  is 
probably  at  least  as  large  an  unknown  factor  to  be  determined 
in  the  equivalent  calculation  still  to  be  done  with  a  kinetic 
physical  model.   The  basic  question,  when  can  transport  co- 
efficients be  calculated  as  averages  over  flux  surfaces,  is 
examined. 
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1.    Introduction 

We  consider  a  macroscopic  resistive  model  subject  to 
Ohm's  law  in  the  form 

E  +u  X  B  =  riJ  .  (!•!) 

The  value  of  t]   is  left  open;  it  can  be  classical  or  "anomalous" 
(and  some  predictions  turn  out  to  be  Independent  of  the  value 
of  T][l,2,3]).   To  complete  the  system,  we  can  take  conservation 
of  mass,  momentum,  and  flixx  in  the  form 

|£+div(pu)  =  0  (1.2) 

7p  =  J  X  B  (1.3) 

SB 


7t 


+  curl  E  -  0  .  (1-^) 


The  isothermal  assumption  and  absence  of  mass  source  implicit 
in  (1.2),  also  the  scalar  pressure  in  (1.3)  and  scalar 
resistivity  in  (1.1)  are  easily  generalized.   But,  in  order 
to  decouple  diffusion  from  MHD  wave  motion,  the  inertia  term, 
pdu/dt,  missing  from  (1.3),  must  remain  excluded. 

This  is  a  much  simpler  physical  model  than  the  ones  treated 
in  most  of  the  recent  literature  which  couple  energy  flow  with 
plasma  diffusion,  usually  in  a  kinetic  guiding  center  or  drift 
model  (e.g.  ['^,5]).   But  it  was  only  recently  recognized  [1] 
that  the  standard  calculation  of  diffusion  in  even  the  simplest 
macroscopic  version  (1.1)  is  woefully  inadequate.   It  is  based 
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on  a  simplification,  viz.  curl  E  =  0  (the  special  theory),  which 
turns  out  to  eliminate  from  consideration  all  but  a  very  special 
class  of  plasma  equilibria.   It  can  be  shown  [1,2,3]  that  this 
special  class  of  equilibria  sometimes  takes  over  after  a  suffi- 
cient lapse  of  time  (probably  too  long  to  be  relevant  in  present 
experimental  Tokomaks),  but  more  often  the  special  theory  is  not 
representative  of  the  general  case.   Independent  of  the  quanti- 
tative relevance  of  the  macroscopic  model  (1.1)  to  a  hot  plasma 
with  large  mean  free  path,  its  qualitative  understanding  is 
crucial  to  the  understanding  of  kinetic  models  for  which  the 
theory  is  still  at  the  earlier  stage  of  the  standard  macro- 
scopic calculation  with  curl  E  =  0.   The  neoclassical  orbit 
considerations  supply  some  of  the  physics  which  is  missing  in 
the  macroscopic  model;  but  the  additional  physics  which  is  miss- 
ing because  of  the  specialized  profiles  in  the  kinetic  model  is 
likely  to  yield  equally  important  quantitative  corrections. 

The  basic  reason  for  the  complexity  of  classical  diffusion 
(also  energy  flow)  is  that  what  is  called  a  "diffusion  (heat 
transfer)  coefficient"  is  not  a  traditional  transport  coefficient. 
It  is  not  a  local  property  of  the  plasma;  rather  it  is  a  "black 
box"  solution  of  a  global  boundary  value  problem.   We  can  com- 
pare the  elementary  concept  of  a  local  viscosity  coefficient 
with  the  drag  coefficient  for  an  airfoil,  involving  global 
interactions  of  wakes,  boundary  layers,  shock  waves,  etc.  — 
each  one  dependent  on  viscosity;  plasma  diffusion  is  more 
akin  to  the  drag  coefficient  than  to  the  viscosity  coefficient. 
The  special  theory  (curl  E  =  0)  can  frequently  be  reduced  to 
a  calculation  of  local  transport  coefficients;  the  general 
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theory  almost  never  can. 

This  distinction  is  not  mere  semantics.   For  example, 
all  practical  numerical  calculations  of  Tokomak  behavior  are 
done  in  a  straight,  cylindrical  geometry  with  a  diffusion  (or 
heat  transfer)  coefficient  which  (among  other  ad  hoc  factors) 
is  theoretically  corrected  for  the  toroidal  curvature  (e.g. 
[6,7])-   But  the  correction  factor  is  taken  from  a  global, 
geometry-sensitive,  convective  calculation  (such  as  [4]  or 
[5]);  moreover  the  theoretical  calculation  is  valid  only 
after  complete  relaxation  of  a  transient  —  and  it  is  the 
transient  which  is  the  object  of  study  of  the  numerical 
calculation.' 

Without  inertia  in  (l.3)j  at  least  two  time  scales 
appear,  separated  by  three  or  more  orders  of  magnitude.   With 
inertia,  faster  MHD  time  scales  appear,  increasing  the  spread 
by  at  least  two  more  orders  of  magnitude.   No  numerical  cal- 
culation which  includes  the  fastest  scale  (MHD  waves  and 
compression)  can  include  the  slowest.   But  elimination  of 
inertia  makes  the  system  very  nonstandard  with  constraints 
replacing  time  derivatives;  thus  elaborate  theoretical  pre- 
paration is  required  before  one  can  proceed  with  a  numerical 
calculation  [3]. 

Basic  features  of  the  general  theory  were  described  in 
[1];  some  specific  Tokomak  applications  have  been  given  in 
[2];  details  of  the  general  theory  will  appear  [3]-   In  this 
paper  we  briefly  summarize  these  results  and  discuss  the  over- 
all significance  of  constraints  and  averaging  on  transport 
properties . 
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2.    Critique  of  the  Special  Calculation 

Almost  all  diffusion  calculations  (both  macroscopic  and 
kinetic)  are  "based  on  two  postulates: 

(i)   small  resistivity,  implying  slow  passage  through  a 
family  of  quasistatic  equilibria  (Vp  =  J  x  B  or  an  equivalent 
guiding  center  equilibrium) . 

(ii)  low  P,  implying  curl  E  =  0. 
Although  these  seem  eminently  plausible,  after  examination 
they  turn  out  to  be  unduly  restrictive. 

There  are  several  consequences  of  these  assumptions,  of 
which  we  list  three: 

(a)  the  net  flow  can  be  calculated  explicitly  and  very 
easily 

(b)  the  pertinent  equilibria  are  restricted 

(c)  the  conclusions  are  very  sensitive  to  slight  changes 
in  the  model. 

Implicit  in  postulate  (i)  above  is  that  there  exist s_  a 
family  of  neighboring  equilibria  (through  which  the  plasma 
slowly  diffuses).   This  is  not  necessarily  so  because  of  the 
phenomenon  of  marginal  equilibrium  [8].   For  example,  a  mirror 
machine  which  is  bent  into  a  bumpy  torus  may  require  an  abrupt 
change  in  the  topology  of  the  pressure  surfaces  (this  phenomenon 
is  distinct  from  the  formation  of  magnetic  islands  by  loss  of 
symmetry) .   Similarly,  in  a  straight  circular  screw  pinch, 
certain  helical  perturbations  nius_t  change  the  pressure  topo- 
logy.  In  addition  to  its  effect  on  diffusion,  this  also  re- 
quires reexamination  of  stability  concepts  such  as  resistive 
tearing  "modes".   Standard  estimates  of  resistive  growth  rates 


-4- 


are  Irrelevant  since  they  do  not  describe  an  Instability  at 
all  but  a  change  In  equilibrium  topology;  to  follow  the 
consequences  of  the  change  in  topology  requires  a  nonlinear 
calculation  of  an  initial  value  problem  rather  than  normal 
modes  or  growth  rates. 

With  regard  to  point  (a),  the  calculation  is  miraculously 
able  to  ignore  all  the  fine  details  of  the  global  interplay 
between  convection  and  diffusion,  and  it  can  always  be  carried 
out,  both  when  it  is  correct  and  incorrect.   For  example,  it 
is  known  that  neither  scalar  pressure  nor  guiding  center 
equilibria  exist  in  asymmetric  geometries  [9];  but  the  calcu- 
lation, nevertheless,  yields  a  diffusion  rate  [10,11];  (this 
is  not  an  abstract  objection  —  it  is  easily  verified  that  the 
"solutions"  do  not  satisfy  the  equations). 

With  regard  to  consequence  (b)  above,  when  there  is  no 
driving  emf  the  mean  parallel  current  on  each  magnetic  surface 
must  vanish  [12] 

<riJ.B>  =  0  (2.1) 

More  generally  [3]^  in  a  Tokomak  or  multipole  geometry  with 
driving  emf's,  the  restriction 

<TiJ.B>-  c^^^'(V)  +C2V'2'(V)  (2.2) 

must  be  satisfied  where 


c  ^6     E-dx  (2.3) 
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are  the  two  (poloidal  and  toroidal)  constant,  path-independent, 
emf's,  and  ''P-,    =   ip   and  ''P^^   are  the  poloidal  and  toroidal  flixxes. 

The  general  static  equilibrium  is  characterized  by  two 
arbitrary  profiles,  e.g.  p(V)  and  V'(V),  or  p(^)  and  q(^)  =  ^2'(^) 
(q  is  the  reciprocal  rotation  number  or  Tokoraak  "stability" 
parameter)  .   The  standard  calculation,  based  on  assux'iptions 
(i)  and  (ii),  is  automatically  limited  to  a  single  arbitrary 
profile  through  the  conditions  (2.1)  or  (2.2).   A  natural 
question  is  whether  these  special  profiles  are  representative 
of  the  general  case.   The  answer  is  that  the  evolution  in  time 
(including  diffusion  properties )ls  even  qual i t at i ve ly  different 
for  the  special  and  general  profiles,  both  physically  and 
analytically.    The  analytical  features  have  necessitated 
development  of  an  elaborate  and  entirely  nonstandard  mathematical 
structure  [3]-   As  just  one  example,  the  classical  boundary 
condition  that  at  a  fluid  interface  the  normal  component  of 
fluid  velocity  equals  the  velocity  of  the  interface  is  violated 
in  this  theory. 

Some  of  the  physical  conclusions  arising  in  the  general 
theory  which  may  be  more  immediately  relevant  to  Tokomak 
experiments  have  been  s\jmmarized  in  Ref.  [2]: 

(1)  present  Tokomaks  are  inherently  transient  since 
they  do  not  persist  long  enough  to  settle  down;  this  casts 
some  doubt  on  the  significance  of  scaling  laws. 

(2)  the  fact  that  the  plasma  does  not  appear  to  be 
transient  is  probably  due  to  the  dominant  effect  of  the  neutral 
influx. 

(3)  noncircular  plasma  cross-section  can  give  a  decided 


advantage  with  regard  to  diffusion  rate  and  P  (a  stability 
advantage  has  "been  noted  previously  [13]). 

This  transient  conclusion  was  seemingly  disputed  in  a 
later  calculation  [1'-+].   There  is  actually  no  contradiction 
[15] J  since  a  steady  state  is  a^ssumed  in  Ref.  l4  which  is  an 
ad  hoc  calculation  based  on  several  experimental  inputs  which 
would  undoubtedly  change  if  the  experiment  were  extended 
sufficiently  in  time. 

With  regard  to  point  (c),  the  sensitivity  of  the  standard 

calculation,  we  point  to  at  least  three  distinct  sources.   The 

classical  skin  effect  has  a  diffusion  coefficient  D   ~  ri/iJ,  . 

o    '   o 

The  "classical"  plasma  diffusion  coefficient  is  approximately 

D-,  ~  ^D  .   Diffusion  in  a  general  equilibrium  is  characterized 

-3 
by  a  competition  between  these  two  processes.   Since  p  ~  10 

_  o 

in  a  Tokoraak  and  possibly  10   in  an  alkali  plasma,  even  a 
very  slight  coupling  with  the  much  faster  skin  effect  might 
raise  D-,  enormously;  (it  is  surprising  that  only  factors  of 
one  or  two  orders  of  magnitude  have  so  far  arisen  in  the  theory) . 

Another  cause  of  great  sensitivity  is  the  complexity  of 
the  convective  flow.   Comparing  the  tx\j'o  different  weighted 
averages 

<u-Vv>>  ~  net  flow 

(2.4) 
<r^><u-V^/r^>  , 

we  should  expect  them,  to  be  approximately  equal  in  a  large 
aspect  torus  ( 6r/r  is  small  on  a  flux  surface),  but  instead 
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o 
we  find  that  their  ratio  is  a  large  number  (~q'~  where  2.5  <  q  <  ^1- 

in  a  Tokomak) .   In  a  finite  aspect  torus,  where  the  two  expressions 

can  be  expected  to  be  somewhat  different,  they  are  even  found 

to  scale  differently.   The  point  is  that  convection  is  strongly 

geometry  dependent  (whereas  a  true  transport  coefficient  would 

not  be) . 

Related  to  this  phenomenon  is  that  fact  that  D,  for  a 

2 
straight  circular  cylinder  is  increased  by  a  factor  1  +  2q 

(Pfirsch-Schluter  [l6])in  a  torus  of  arbitrarily  large  aspect 

ratio  even  though  q  itself  is  given  the  same  value  in  the 

straight  and  toroidal  geometries. 

A  similar  calculation  with  resistivity  ri  ^  r\{i/)    varying 
slightly  (i.e.  of  order  l/aspect)  on  a  flux  surface  gives  a 
correction  factor  to  D,  which  may  be  large  compared  to  unity 
[3,15]  just  as  for  the  Pfirsch-Schluter  factor. 

Another  independent  cancellation  occurs  in  a  Tokomak  (or 
in  any  discharge  with  6E'dx  j^  0) .      In  a  passive  plasma,  the 
net  diffusion  is  directly  related  to  the  dissipation  <tiJ  >; 
in  a  discharge  it  is  smaller  by  an  order  of  magnitude.   A 
similar  sensitivity  to  the  presence  of  a  driving  emf  has  been 
found  for  large  mean  free  path  in  the  form  of  trapped  particle 
pinching  [17]. 
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3.    The  General  Classical  Calculation 

With  a  general  static  equilibrium  profile,  the  assumption 

curl  E  =  0  is  not  justified  even  with  arbitrarily  small 

resistivity  and  P.   One  contributive  factor  is  that  the  skin 

effect  (field  penetrating  plasma)  and  plasma  diffusion  (plasma 

penetrating  field)  are  coupled  through  the  pressure  balance, 

giving  rise  to  a  complex  convective  flow  which  modifies  both 

primitive  concepts.   The  skin  effect  in  a  deformable  plasma 

is  not  the  same  as  in  a  resistive  solid  conductor,  since  this 

would  clearly  violate  the  pressure  balance.   The  plasma 

diffusion  which  accompanies  the  skin  penetration  is  entirely 

nonstandard,  even  in  scaling.   There  is  a  very  strong  geometrical 

dependence;  circular  sections  are  qualitatively  different  from 

noncircular;  Tokomaks  are  qualitatively  different  from  hard 

cores  or  multipoles;  there  are  singular  limits  as  the  toroidal 

or  poloidal  field  dominates  the  other;  the  skin  effect  can 

involve  one  field  component  or  both,  and  in  a  9-pinch  or  bumpy 

torus  there  is  no  skin  effect  at  all  on  the  scale  D  .   Not 

o 

only  the  physical  conclusions,  but  the  mathematical  structures 
differ  in  these  cases.   Entirely  different  numerical  schemes 
are  required  in  order  to  "march"  in  different  cases,  and  a 
physical  variable  (such  as  f)    which  is  advanced  because  of  a 
time  derivative  in  one  case  (^f/dt)    will  advance  in  virtue  of 
an  implicit  constraint  in  another  case. 

In  a  large  number  of  cases  a  standard  s e que nc e ,  skin 
effect  followed  by  "classical"  diffusion  can  be  verified 
theoretically.   For  example,  in  a  low  P,  large  aspect  Tokomak, 
the  toroidal  current  will  decay  quickly  in  a  time  related  to 
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D  .   This  process  v;ill  approximate  the  standard  skin  effect 
if  the  plasma  cross-section  is  approximately  circular;  it  will 
have  the  same  scaling  but  will  differ  quantitatively  from  the 
standard  skin  effect  if  noncircular.   During  this  process  the 
pressure  contours  follow  the  changes  in  flux  surface  shape; 
the  pressure  remains  approximately  constant  on  a  contour  of 
fixed  volume.   The  toroidal  field  cannot  be  said  to  diffuse. 
Instead  it  fills  in  the  gaps  to  maintain  pressure  balance. 
During  the  skin  diffusion  the  plasma  does  diffuse  slightly 
relative  to  a  fixed  volume,  but  at  a  nonstandard  rate.   The 
diffusion  rate  is  not  local  and  the  process  cannot  be  described 
in  terms  of  mean  values.   After  a  number  of  skin  times,  the 
system  settles  down  to  one  of  the  special  profiles  of  the 
standard  calculation.   Provided  that  certain  additional  con- 
ditions are  met,  the  system  then  diffuses  slowly  to  a  imiversal 
equilibrium  with  no  adjustable  parameters,  independent  of  the 
initial  profile  (and  even,  to  some  extent,  independent  of  the 
value  of  T]    [2])  . 

This  standard  sequence  is  not  necessarily  followed,  e.g. 
at  finite  P  or  even  at  low  ^  at  finite  aspect.   The  ultimate 
equilibrium  may  not  be  reached  because  of  instability  (this 
is  a  new  type  of  resistive  instability),  or  it  may  not  be 
reached  because  it  does  not  exist.   For  example,  a  Tokomak 
equilibrium  in  a  straight  circular  cylinder  can  be  found  with 
a  suitable  mass  source  (ionization  of  neutrals)  taking  the 
place  of  the  discharge,  i.e.  with  AE.dx  =  0  [2].   But  no  such 
equilibrium  exists  for  arbitrarily  small  toroidal  curvature  [3]. 
This  is  significant  with  regard  to  the  "bootstrap  Tokomak" 
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concept  [18,19]. 

By  varying  the  plasma  cross-section,  the  rate  of 
approach  to  the  universal  equilibrium  can  "be  accelerated, 
and.  the  plasma  can  be  given  more  interesting  parameters 
(e.g.  higher  P)  [2]. 

A  beginning  of  a  transient  (i.e.  of  a  general)  theory 
can  be  found  in  the  guiding  center  model  as  the  "trapped 
particle  pinch  effect"  [17]- 

The  details  of  the  diffusion  theory  for  a  general 
profile  are  quite  elaborate  [3];  one  qualitative  reason  for 
the  vast  difference  between  the  special  and  the  general  theory 
can  be  described  in  terms  of  their  relation  to  irreversible 
thermodynamics,  which  we  now  discuss. 
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4.    Transport  Coefficients,  Symmetry  and  Averaging 

Perhaps  one  reason  that  the  general  theory  of  diffusion 
for  the  classical  system  (1.1)-(1.4)  remained  undeveloped  for 
so  long  is  that  it  cannot  be  put  into  a  standard  "transport 
coefficient"  framework  as  caxi  the  special  theory.  It  turns 
out,  however,  that  even  the  special  theory  becomes  standard 
only  with  further  restrictions.  This  is  best  seen  after  we 
discuss  several  simpler  examples. 

All  standard  transport  coefficients  are  subject  to  the 
Onsager  symmetry  relations.   However,  there  are  additional 
nonthermodynamic  symmetries  which  arise  from  quite  a  different 
source,  as  for  example  through  a  Green's  function  or  topological 
symmetry.   Entirely  new  problems  arise  when  forces  or  fluxes 
are  made  subject  to  constraints  or  are  averaged.   Also,  the 
extension  of  the  symmetry  principle  to  include  skew  symmetry 
for  terms  linear  in  B  (Hall  effect)  may  not  hold  for  global 
or  averaged  quantities. 

Consider  first  the  classical  problem  of  the  flow  of  heat 
in  a  simply  connected  solid  body  (Fig.  1) .   We  single  out 
portions  of  the  boundary,  S-,,S2,  ...  which  are  held  at  temper- 
atures T, ,Tp, ...   Solutions  of  the  elliptic  system 

q  =  -AVT,    div  q  =  0  (4.1) 

can  be  described  by  a  symmetric  Green's  function.   Thus  the 
net  heat  flow 

\=  f    q-dS  (4.2) 

^i 
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Is  related  to  the  set  of  temperatures  T.  tiy  a  symmetric  heat 
conductivity  matrix 


Q.  =  A.  .T.  (4.3) 


The  original  problem  has  a  single  local  transport  coefficient 
A  and  consequently  has  no  Onsager  identities.   The  discrete 
formulation  in  terms  of  T.  and  Q.  gives  an  array  of  forces 
and  fluxes,  and  the  sjrrametric  relation  (4.3)  arises  from  the 
symmetry  of  the  Green's  function. 

Consider  next  the  dissipative  relation 

E  =  TiJ  (4.4) 

in  a  hollow  toroidal  copper  shell  (Fig.  2) .   We  suppose  that 
ri(x)  is  given  and  that 

curl  E  =  0,    div  J  =  0  (^-5) 

and  J   =  0  on  the  boundary.   This  is  a  classical  elliptic 
system,  and  the  manifold  of  solutions  is  known  to  be  a  two- 
dimensional  linear  vector  space  [20].   A  unique  solution  is 
determined  by  specifying  a  pair  of  periods,  either  emf's 


E«dx  =  c.  (^-6) 

^i 


or  currents 
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/. 


J.dS  =  I.  (4.7) 

Si 


(see  Fig.  2),  which  are  related  by  the  symmetric  resistance 
matrix. 


c.  =  R.  .1.  (4.8) 

1    ij  J 


As  in  the  previous  example,  a  single  local  transport  coefficient 
ri(x)  proliferates  to  an  array  of  forces  and  fluxes  when  examined 
globally. 

In  an  abstract  formulation. 


(4.9) 


where  v.  are  the  "fluxes",  u.  the  "forces",  and  A   is  any 

nonsingular  matrix,  we  can  symmetrize  trivially  by  choosing 

appropriate  new  fluxes  v. .  Given  any  symmetric  nonsingular 
L.  .,  we  simply  set 


V.  =  L.  .u.  (4.10) 

1    ij  J 


where  the  new  fluxes  v  are  defined  in  terms  of  the  old  by 
V  =  la"  V  (there  are  only  minor  modifications  necessary  if 
the  forces  and  fluxes  are  not  independent  and  A  and  L  are 
singular).   Thus  one  can  always  symmetrize.   The  problem  is 
whether  this  occurs  through  a  natural  or  evident  choice  of 
variables  in  a  given  physical  problem. 

The  proper  abstract  formulation  of  Irreversible  Thermo- 
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dynamics  is  to  take  forces  and  fluxes  as  dual  linear  vector 
spaces  (in  the  standard  terminology  of  linear  algebra)  with 
the  entropy  production  as  the  defining  inner  product.   The 
question,  whether  an  abstract  inner  product  (u, v)  is  symmetric 
has  no  meaning.   The  symmetry  of  a  specific  matrix  representa- 
tion of  (u, v)  depends  on  the  choice  of  base  vectors  in  u  and 
in  V.      Similarly,  we  have  just  seen  that  an  abstract  mapping 
from  u  to  V,  such  as  represented  by  (^.9)  or  (4.10),  can  be 
symmetric  or  not  depending  on  the  choice  of  base  vectors. 
However,  the  syimnetry  of  a  map  from  u  to  v  i^elative  to  a  given 
inner  product  (u, v)  does  have  an  invariant  meaning.   This  is 
the  content  of  the  Onsager  principle,  that  the  relation  between 
forces  and  fluxes  is  symmetric  relative  to  a  specific  bilinear 
form,  the  entropy  production. 

If  in  the  examples  above  the  medium  is  anisotropic  and 
there  is  a  matrix  relating  the  three  components  of  J  and  E, 
or  if  there  are  several  coupled  vector  fluxes  J  and  forces 
E  ,  then  we  may  have  local  Onsager  identities  at  a  given  point 
X  of  the  medium.   The  global  symmetry,  relating  emf  to  total 
current,  is  a  consequence  of  two  distinct  symmetries,  one 
thermodynamic  and  local,  and  the  other  primarily  topological. 

The  geometrical  syinmetry  can  be  separated  out  and  is  also 
related  to  an  inner  product.   In  a  domain  of  arbitrary  topology, 
the  linear  vector  space  of  closed  curves  C.  is  isomorphic  to 
the  linear  vector  space  of  cross  cuts  3.  (Alexander  duality 
theoremi,cf  [20]).   A  reciprocal  basis  can  be  found  so  that 


(C.,S  .)  =  5.  .  (^.11) 
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where  (C.,S.)  is  the  algebraic  number  of  Intersections  of  the 
curve  C.  with  the  surface  S.. 

Given  two  vector  fields  X  and  Y  satisfying 


curl  X  =  0 

(4.12) 
div  Y  =  0 


in  the  domain  and  Y  =  0  on  the  boundary,  an  easy  calculation 
gives 


n 


(X,Y)  =  /x-Y  dV  =  Y1   ^.-^^  (^-13) 

^  1 


where 


^.  =  (f  X-dx,    n^  =  r  Y-dS  (^-1^) 

(many  generalizations  and  extensions  of  this  formula,  Hodges ' 
Theorem,  are  given  in  [20]).   The  system  (4.12j  is  incomplete 
and  does  not  determine  X  or  Y.   To  complete  the  system  v;e  need 
a  linear  relation  between  X  and  Y.   The  simplest,  X  =  Y,  deter- 
mines a  manifold  of  solutions  which  is  an  n-dimensional  linear 
vector  space  spanned  by  either  |.  or  rj  .   If  Z  is  a  specific 
solution,  its  dual  representations 


Z  =  E^^x^  =  Et^^Y.  (^.15) 


where  X   and  Y.  are  reciprocal  bases  relative  to  the  inner 
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product    (4.13), 


(x\y^.)    =   &j  (4.16) 


gives    the   syiniaetric   matrix 


e^    =    (Y.,Y    )ti'^'  (4.17) 


Talcing  i.    =  c.  and  n  .  ==  1.,    we  have  proved   that  R.  .  -=  (Y.,Y.) 
is  symtaetric. 

Any  local  symiietric  linear  relation  between  X  and  Y 
which  is  used  to  complete  the  system  will  combine  w'.th  the 
topological  symiietry  to  ^ive  a  globally  syni'iietric  resalt. 

In  a  simple  domain,  dropping  the  boundary  condition  Y  =  0 


and  setting  X  =  V(t),  we  have 


(X,Y)  =  -J<I>Y^   dS  (4.18) 

This  symmetric  inner  product  between  boundary  values  of  *  and 
Y   [which  is  the  source  of  the  symmetry  (4.3)],  coupled  with 
a  symmetric  local  (Onsager)  relation  between  Y  and  VO  gives 
a  globally  symmetric  structure. 

To  illustrate  the  effect  of  averaging,  consider  an  abstract 
mathematical  model  of  forces  u. (x)  and  fluxes  v. (x)  related  by 
a  symmetric  matrix  L. .(x) 

v^(x)  =  L  .(x)  u  (x)  (4.19) 
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If  u .  =  const,  (independent  of  x) ,    then  averaging  clearly 
preserves  the  symmetric  structure, 

jTv^dx  =  {J-L^.6x)u^  (4.20) 

The  same  is  true  if  v.  =  const,  in  terms  of  the  reciprocal 

matrix,  L~  .   The  conditions  u.  =  const,  or  v^  =  const,  are 

natural  and  not  arbitrary.   For  example,  minimizing  the 
quadratic  form 

(u,v)  -/v^l5;J  v.dx  (4.21) 

subject  to  fixed  total  flux 

Tv.dx  -  V.  (4.22) 

gives,  by  the  Lagrange  multiplier  rule, 

seTv-L"-":  v.dx-A.  Tv  dx]  =  0  (4.23) 

the  natural  minimization  condition 


2u.  =  A.  =  const.  (4.24) 

11 


It  is  easy  to  see  that  there  is  no  natural  averaging  which 
preserves  the  linear  structure  when  both  forces  and  fluxes 
depend  on  x. 

We  can  now  understand  the  successful  averaging  in  the 
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first  two  examples  given  above.   The  heat  flow  problem  ('^.1) 
gives  a  simple  symmetric  averaged,  structure  ("^.3)  when  the 
T.  are  taken  to  be  constant  on  each  boundary  segment;  other- 
wise it  does  not.   Similarly,  the  electrical  conduction 
problem  (4.4),  (4.5)  gives  a  simple  averaged  result  because 
curl  E  =  0  guarantees  that  the  forces  c.  are  constant  (path- 
independent)  . 

Averaging  occurs  in  a  more  complex  way  if  J  is  contrained. 
to  lie  in  given  surfaces,  or  if,  because  of  highly  anisotropic 
heat  conductivity,  T  is  constrained  to  be  constant  on  specified 
surfaces.   The  standard  transport  formulation  must  be  modified 
when  there  is  a  constraint  and,  as  we  might  expect,  the  inner 
product  gives  the  clue. 

If  T  is  constrained  to  be  constant  on  a  given  set  of 
surfaces,  T  =  T(-A),  the  inner  product  is 

/vT-q  dV  =  c&  q.dS  -  jT   div  q  dV  (^-25) 

With  T  =  T{f) ,    the  volume  term  vanishes  if  <(div  q/'^-^  0  on 
each  ^-surface.   Thus  restricting  the  admissible  T  forces  us 
to  expand  the  admissible  q  from  div  q  ^=  0  at  each  point  to 
\  div  q^  =  0  on  each  surface. 

If  J  is  contrained  to  lie  in  a  specified  surface, 
V^'J  =  0,  we  can  write 

J  =  V^  X  VcD  (^-26) 

and  evaluate 
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jE'J   dV  =  -/E.curl(cuV^)dV 

(4.27) 
=  -/a)(V^.curl  E)  dV  +  boundary  terms 

Thus  the  restriction  on  J  requires  us  to  relax  E  from  curl  E  =  0 
to  V^.  curl  E  =  0;  E  is  a  surface  gradient  on  each  lA-surface. 
Note  that  taking  the  conventional  equations  div  q  =  0  or  curl  E  =  0 
is  incorrect  (and  overdetermines  the  system)  when  there  is  a 
constraint  such  as  T  =  T(^)  or  V^.J  =0.   (In  the  plasma  dif- 
fusion problem  we  are  permitted  to  retain  curl  E  =  0  rather 
than  V^.curl  E  =  0  only  because  Oh-n's  law  is  relaxed  to  include 
an  adjustable  term  u  x  B.) 

The  constrained  current  problem  is  formulated  as 

E  =  TiJ,    V^.  curl  E  =  0 

(4.28) 
div  J  =  0,    V^.J  =  0 

On  each  surface  we  have  an  elliptic  system  and  (4.8)  generalizes 
to 

c.(^)  =  R.  (-^)I  •(^)  (4.29) 


where  I.(^)  is  the  total  current  within  the  ^-surface.   A 

J 

simple  symmetric  structure  obtains  on  each  ^-surface  since 
the  forces  c.  are  constant  on  a  surface.   But  no  such  structure 
exists  when  averaged  over  ^  unless  curl  E  =  0  rather  than 
V^.curl  E  =  0. 
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Up  to  now  we  have  taken  a  restricted  Onsager  principle 
without  a  magnetic  field.   With  Hall  effect,  the  relation 
between  forces  and  fluxes  is  no  longer  symmetric.   It  is  con- 
ventional to  include  the  skew  symmetry  of  terms  linear  in  B 
(this  is  not  a  physical  principle  at  all,  but  a  consequence 
of  the  fact  that  B  is  not  a  vector)  as  an  extension  of  the 
"symmetry"  principle.   This  convention  is  not  tenable  when 
we  turn  to  global  considerations.   A  Kail  term  in  Olim's  law, 
E  =  qj  +  aj  X  B  or  in  heat  conduction,  q  =  -AVT  +  xq  x  B,  makes 
the  differential  system  non-self  adjoint,  destroys  the  symmetry 
of  the  Green's  function,  and  usually  destroys  the  global 
symmetry  between  forces  and  fluxes.   Symmetry  cannot  be  re- 
stored by  any  expedient  related  to  reversal  of  part  of  the 
Uiatrix,  since,  for  example,  various  magnitudes  and  orienta- 
tions of  B  may  be  averaged  globally. 

The  simplest  example  with  Hall  effect  is  two-dimensional 
heat  flow;  we  take  the  same  problem  as  before.  Fig.  1.   We 
have 

q  =  -AVT  +Tq  x  B  (^.30) 

where  B  is  a  constant  vector  perpendicular  to  the  plane  of 
interest,  and  A  and  x  are  constant.   Using  div  q  =  0,  we  can 
set 

q  =  n  X  V^  (^.31) 

where  n  is  the  transverse  unit  vector.   We  have 
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n  X  V^  =  V*  (^.32) 


where 


t)  =  tB^A-  at  (^.33) 

From  (4.32),  we  see  that  <t>   and  f   are  conjugate  harmonic 
potentials  related  by  a  symmetric  Green's  function.   In  par- 
ticular, *  =  const,  on  each  segment  S.  of  the  boundary  gives 
a  symmetric  global  relation  between  <i>.    and  Q.  .   But  <t)  =  const, 
on  S.  is  inconsistent  with  T  =  const,  {f   must  vary  if  there 
is  any  heat  flow  through  S.).   In  fact,  we  have  symmetrized 
by  changing  bases  in  (4.33) (as  we  have  shown,  this  can 
always  be  done).   Even  without  making  the  problem  discrete, 
i.e.  without  averaging,  we  can  verify  that  the  boundary 
values  q  (S)  are  not  symmetric  in  terms  of  T(S),  although 
they  are  symmetric  in  terms  of  '''(S). 

It  should  be  pointed  out  that  when  a  Hall  term  sach  as 
q  X  B  points  in  an  innocuous  direction,  the  desymmetrizing 
effect  may  disappear.   For  example,  take 

E  =  iiJ  +aj   X  B  (4.3^) 

together  with  the  remainder  of  (4.28) (note  that  B  is  a  given 
constraining  field  with  B-Vi;//  =  0,  and  in  this  context,  has  no 
relation  to  J)  .   The  surface  equation  is  the  s.-ame  as  before, 
and  the  normal  component  of  (4.34)  merely  serves  to  determine 
the  component  E-V^^  after  J  has  been  found.   Note  that  this 
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system  is  not  self  adjoint,  yet  (4,29)  remains  symmetric. 
More  generally,  a  set  of  fluxes  q  and  forces  VT  [f) 

gives  a  symmetric  averaged  matrix  even  with  Hall  effect. 

ct 
The  same  is  true  of  a  set  of  fluxes  J  each  constrained  by 

a 
Vip.  J     =   0.   However,  a  system  with  mixed  constraints,  coupling 

(q,  T),T  =  T(^)  with  (J,  E),J'Vv^  =  0,  is  not  syimnetric  and  does 

not  allow  averaging  in  the  presence  of  Hall  effect. 

Turning  to  the  plasma  diffusion  problem,  since  there  is 

only  one  local  transport  coefficient,  r\,    any  global  symmetries 

that  arise  are  not  basically  thermodynamic.   With  the  special 

equilibrium  (2.2),  the  dissipation  takes  the  form 

Aj^dV  =  c^I-^  +C2I2-  r(u.Vp)dV  (^-35) 

where  I-,  and  Ip  are  the  total  poloidal  and  toroidal  currents. 
An  unusual,  symmetric  "matrix"  arises  [3]  relating  the  forces 


c^,C2,  -p'(V)  (4.35) 


and  the  fluxes 


I,,I^,  I   u.dS  (4.37) 


-l'"2 


Note  that  c-,,Cp  and  I-i  j  Ip  are  global  (averaged  over  the  entire 
domain),  whereas  p'(V)  and  ou-dS  are  local  to  a  surface.   The 
symmetry  under  averaging  is  preserved  because  c-,  and  Cp  are 
constant  (curl  E  =  0)  and  p  is  constant  on  a  flux  surface. 
Variation  of  11  with  tp   or  even  over  a  flux  surface  is  irrelevant 
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to  the  structure  (but  can  change  the  numerical  values 
considerably) . 

It  is  immediately  evident  that  inclusion  of  heat  flow, 
either  by  itself  or  coupled  to  diffusion,  can  yield  a  standard 
structure  after  averaging  over  fliix  surfaces  only  when  T  =  T(V'). 
This  restriction  is  inadmissible  in  the  present  theory,  but 
it  may  be  valid  at  very  large  mean  free  path.   Coupled  heat 
flow  and  diffusion  yields  an  Innocuous  Hall  effect  (both  are 
of  type  q"  above).   But  in  a  discharge,  6E-dx  ^   0,  symmetry 
is  lost  and  there  is  no  longer  any  simple  relation  between 
averaged  forces  and  fluxes  (since  the  system  now  has  mixed 
constraints) . 

Similarly,  with  fluid  flow  (steady  state  rather  than 
static  equilibrium)  and  Hall  effect,  since  the  nondlssipative 
equilibrium  has  distinct  flux  and  pressure  surfaces,  p  -^  p{f) 
[21],  even  diffusion  by  Itself  cannot  be  described  in  terms 
of  net  flow  and  other  flux  surface  averages.   In  a  guiding 
center  model,  the  equilibrium  can  sometimes  be  written  [22] 

K  x   B  =    {hip/dTp)S/f 

K  =  curl    (aB)  (^.33) 

a   =   1  -  (Sp/dB)/B 

where  p  =  p(^, B)  is  the  parallel  pressure.   K  takes  the  place 
of  J  and  is  constrained  to  lie  in  "/'-surfaces.   In  Ohm's  law. 
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E  +u  X  B  =  ..  .  (^.39) 

maltiplication  "by  K  gives  an  averaged,  term  \E'K)>  v/hich  can  be 
evaluated  in  terms  of  global  emf's,  but 

<u  X  B-K>=  -  <  (u'V'^)ap/S^>  (4.40) 

is  not  simply  interpretable.   Specializing  to 

P  =  Pi(^)P2(B)  ,  (4.41) 

we  have 

<u  X  B-K>-  -p-^'<(u-V^)p2>  (4.42) 

which  indicates  that  the  force  p.,  '  may  be  associated  with  a 
weighted  plasma  flow.   The  net  plasma  flow  occars  in  an 
averaged  calculation  only  to  the  extent  that  B  =  const,  on 
flux  surfaces  (see  also  [15]). 

In  a  kinetic  model,  there  are  no  local  transport  coeffi- 
cients; it  requires  flux  surface  averaging  to  introduce  this 
concept  at  all.   There  is  the  same  distinction  between  a 
general  and  a  special  theory  (the  latter  with  curl  E  =  0  and 
restricted  profiles).   As  in  the  small  mean  free  path  theory, 
there  are  no  averaged  transport  coefficients  in  the  general 
theory,  and  there  are  strict  limitations  on  the  use  of  averaging 
in  the  special  theory.   The  consequence  is  that  there  is  no 
widely  applicable  theory  of  transport  coefficients,  either 
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locally  or  averaged  over  surfaces.   In  other  words,  each 
problem  is  solved  globally  as  a  special  boundary  value 
problem  for  a  special  profile,  and  there  is  no  direct  carry- 
over to  another  problem  with  different  shape  or  profile  (as 
there  would  be  for  a  legitimate  transport  coefficient) . 

To  elucidate  the  general  theory,  consider  the  relation 
between  conservation  equations  and  dissipation.   Conservation 
of  energy  in  a  static  fluid  can  be  v/ritten 

C^ll  +  divq.^0  (^.^3) 

which  together  with  the  transport  relation 

q  -  -AVT  {^.^'^) 

gives  a  simple  dissipative  diffusion  equation  for  T.   Note  that 
the  force  is  a  gradient  of  the  quantity  which  is  conserved. 
The  isothermal  conservation  of  mass 

|£  +  div  (pu)  =  0  (^.45) 

together  with  "classical"  plasma  diffusion, 

u  =  -D^Vp/p  (4.46) 

gives  a  similar  diffusion  equation  for  p.   Conservation  of  flux. 


ll  +  curl  E  =  0  (4.47) 
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together  with  the  transport  relation  E  =  tiJ,  gives  a  similar 
diffusion  equation  for  flu:<:. 


=  -curl[(riA  )curl  B]  .  (4.48) 


o 


There  is  a  complete  analogy  of  flu^c,  mass,  and  energy  conser- 
vation.  But,  curiously,  the  analogy  holds  only  If  E  is  called 
the  flux  and  J  the  force;  also  the  heat  conductivity  A  takes 
the  same  position  as  electrical  resistivity  t],    viz.  as  a  diffu- 
sivlty  coefficient.   This  has  a  nimber  of  interesting  con- 
sequences which  we  shall  not  pursue  here. 

In  all  these  time  dependent  formulations  the  transport 
relation  itself  does  not  depend  explicitly  on  time.   This  is 
the  crucial  difference,  in  the  present  context,  between  the 
special  plasma  diffusion  theory  and  the  general  theory.   In 
the  latter,  although  the  local  transport  relations  (when  they 
exist)  are  always  time-independent,  the  averages  over  flux 
surfaces  can  never  be  separated  from  the  global  transient 
solution  of  bhe  entire  system.   There  is  no  doubt  that  the 
situation  will  be  found  to  be  the  same  in  all  finite  or  large 
mean  free  path  models  when  the  general  time-dependent  case  is 
attempted.   We  conclude  that  the  special  solutions,  with 
curl  E  =  0,  give  limited  information  about  plasma  transport 
behavior  in  general.   And  even  within  the  special  theory,  only 
in  certain  well-defined  special  cases  do  the  transport  properties 
reduce  to  a  calculation  of  fluix  surface  averages. 
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5.    Conclusion 

In  a  small  mean  free  path  plasma,  local  transport 
coefficients  exist  at  each  point  of  the  plasma.   Transport 
coefficients  in  a  finite  mean  free  path  system  can  only 
exist  as  a  consequence  of  averaging  over  fliix  surfaces.   The 
condition  curl  E  =  0  does  not  follow  from  low  P  "but  from  a 
specialization  of  profiles  (and  sometimes  from  aging  the 
system) .   In  the  special  theory,  with  curl  E  =  0,  in  addition 
to  the  fact  that  this  postulate  represents  only  a  small  part 
of  the  physical  world,  averaging  leads  to  a  conventional 
system  of  transport  coefficients  only  after  additional 
specialization.   Even  when  there  do  exist  averaged  transport 
coefficients,  they  frequently  cannot  be  found  directly  from 
an  entropy  production  quadratic  form  because  of  the  inter- 
ference of  Hall  effect  \\rhich  may  destroy  the  averaged  symmetry. 

The  general  theory,  with  curl  E  7^  0,  exhibits  many  new 
dissipative  properties  and  requires  new  analytic  techniques. 
It  turns  out  that  numerical  calculations  of  transients  must 
be  done  in  an  actual  toroidal  geometry.   Theoretical  results 
have  thus  far  been  obtained  only  for  the  macroscopic  case, 
and  even  here  mostly  for  diffusion,  uncoupled  to  energy  flow. 
It  has  been  predicted  [1,2]  that  corrections  of  up  to  two 
orders  of  magnitude  should  be  expected  to  arise  when  this 
theory  is  generalized  to  finite  or  large  mean  free  path  problems. 


-28- 


figure  1 
Simply  Connected.  Domain 


figure  2 
Toroidal  Domain 
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